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Abstract 

We consider populations structured by a phenotypic trait and a 
space variable, in a non-homogeneous environment. In the case of sex- 
ual populations, we are able to derive models close to existing mod- 
els in theoretical biology, from a structured population model. We 
then analyze the dynamics of the population using a simplified model, 
where the population either propagates through the whole space or 
it survives but remains confined in a limited range. For asexual pop- 
ulations, we show that the dynamics are simpler. In this case, the 
population cannot remain confined in a limited range, i. e. the popu- 
lation, if it does not get extinct, propagates through the whole space. 



1 Introduction 

In this paper, we are interested in populations that are structured by a con- 
tinuous phenotypic trait »6l and a continuous space variable living 
in a non- homogeneous environment: we will consider a phenotypic trait of 
best adaptation 6{x) that depends on the space variable. This type of popu- 
lation have been studied in [2U [HI [25] , using mostly numerical simulations. 
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This type of models can for instance be used to study the distribution 
of a species along an environmental gradient (such as the north-south gradi- 
ent of temperature in the northern hemisphere): To study the range of the 
species and its local adaptation, one should consider both evolution and spa- 
tial dynamics (see [HI [191 ESI E]). Those models are also useful to study the 
impact of an environmental change (e.g. global warming) on a population 
(see [21 USES]). 

Our work is largely based on [T5] , and related articles [2U [121 E51 EE] • In 
[15] . a partial differential equation model describing the spatial and evolu- 
tionary dynamics of a population has been introduced: 

f d t N(t, x) - A x N(t, x) = (1 - \(Z(t, x) - Bxf - N(t, x)) N(t, x), 
\ d t Z(t,x) -A x Z(t,x) = 2d x (log N(t,x))d x Z(t,x) + A(Bx - Z(t,x)). 

(1) 

where N(t, x) is the population density at the location x, Z(t, x) its mean 
phenotypic trait, and A, B are two parameters. Numerical simulations where 
run for this model, and they showed that depending on A and B, three 
biological scenarios were possible: 

• if B is large (the environment changes rapidly in space), then the pop- 
ulation goes extinct, 

• for intermediate values of B, the population survives, but remains in a 
limited area, 

• if B is small, the population invades the whole space. 

In this paper, we show how ([TJ (indeed, the closely related model (TT5|) ) 
can be derived from a structured population model in the case of sexual 
populations, but isn't appropriate to model asexual populations. We also 
introduce a simplified model that allows us to investigate the dynamics of 
the population. 

In Section 2 , we introduce a structured population model for the evolu- 
tion of sexual populations structured by both a phenotypic trait and a space 
variable. To construct this model, we add a spatial variable to a well esta- 
blished local selection-mutation model (similarly to many kinetic models in 
physics or chemistry, see [28]). This structured population model can also be 
seen as a continuous version of the model introduced in [25J. In this struc- 
tured population model, three parameters appear: A, B, C. We show that 
a model very close to (JTJ) can be obtained as a formal limit of our structured 
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population model when C is large (which means that many generations are 
necessary to obtain a significant growth of the population), see (fT3"j) . We 
explain in Rem. [2] why we couldn't obtain exactly the same model. 

Moreover, we consider a formal limit of ( TT3|) when A is small (which 
corresponds to the values of A considered in [19]). This way, we obtain a 
simpler model on Z(t,x) only: 



In Section 3, we analyse the model ([2]) derived in Section 2. Unfortunately, 
this equation has singularities that are obstacles to have a well-defined prob- 
lem: we show that viscosity solutions exist, but are not unique. Nevertheless, 
this simple model allows us to describe two of the three possible scenarios 
from [19]: invasion fronts, and populations remaining in a limited area. The 
extinction phenomena cannot be observed here because of our assumption 
that A is small. 

In Section 4, we investigate the case of asexual populations. Similarly 
as in Section 2, we introduce a structured population model for asexual 
populations structured by both a phenotypic trait and a space variable. After 
a rescaling, we obtain a structured population model that depends on two 
parameters, A, B. We show that in this case, only two biological scenarios 
are possible: either the population goes extinct, or it spreads to the entire 
space (we were however only able to show this last result in a weak sense, 
see Thm. [S]). This result shows in particular that the model ([T]) from [TS] 
does not apply to asexual populations. 

2 Sexual populations 
2.1 The model 

We start from a classical model describing the evolution of a population 
structured by a phenotypic trait only (see e.g. [8j[T5l[23], and [THEEJEO] f° r 
mathematics properties of this kind of models). Let n(t,v) be the density 
of the population at time t > and phenotypic trait v G R. We assume 
that the fitness depends on the square of distance between v and an optimal 
adaptation trait 6, and is altered by the population size. Then, the fitness 
s[n(t, -)}(v) of an individual of phenotypic trait v living among a resident 





d x Z(t,x) 



(2) 
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population n(t,-) is then given by: 



s[n(t, -))(v) = r max - ^r(v - 9) 2 - J n(t, w) dw, (3) 

where r max > denotes the maximal growth rate of the population, V s 
measures the strength of the selection toward the optimal trait 9, and K is 
the carrying capacity of the environment. 

The fitness is well-defined in the case of asexual populations: it is the rate 
of births of offspring of trait x minus the rate of death. In the case of sexual 
populations, however, the situation is more complicated, since reproduction 
requires two parents, of traits t>* and v'^, that give birth to an offspring of 
trait v, usually different from v* or v'^. This has two consequences (see [7]): 

• We then need to define Q(-, v*, v'*) the distribution function of the trait 
of the offspring. We will analyze the properties of Q in Subsection 12.21 

• We have to define the birth rate and the death rate separately. We will 
assume here that the birth rate is a constant, equal to 7 > r max . 

If moreover we do not differentiate male and female (for instance because 
they have the same distribution), assume that mating is random and uniform 
among the population, and that the number of offspring is proportional to the 
population density (the idea being that the number of births is proportional 
to the number of females), then the evolution of the population structured 
by a phenotypic trait only is described by (see [15]): 



d t n(t,v) 



-(7 - r ma x) ~ -j - ^r(v - Of - / //(/. (r )d(r 



n(t, v) 



, 1 n{t,v,)n(t,v'J , , 

+^ / / ~T^l — VI — Q(v,v*,v,)dv*dv m , (4) 
J n{t, w) dw 




where 7 is the birth rate (7 > r max ), that we assume constant within the 
population (the selection occurs in the death term), and y the additional 
death rate due to lethal mutations (see [12]). 

We consider next populations that are structured by a phenotypic trait v 
as above, but also by a space variable x G K: n(t,x,v). We assume that the 
selection-mutation process described above occurs locally in space, but that 
individuals move randomly in space (which we model by a diffusion of rate 
a x ), and that the trait of optimal adaptation changes linearly in space: 

6{x) = bx. (5) 
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We then get the following model for sexual populations: 

d t n(t, x, v) - alA x n(t, x, v) 



T 1 r 

il-r max )---— (v-bx) - — 



n(t, x, w) dw 



+ 7 




n(t, x, v*)n(t, x, vl] 



f n(t, x, w) dw 
For an existence theory for we refer to 



Q(i>,f*X) dv* dv*. 



n(t, x, v) 
(6) 



2.2 Properties of the sexual reproduction kernel Q 

In Subsection 12.31 we will rescale ([6]) to show that it indeed only depends 
on three parameters. But to be able to do so, we need first to define more 
precisely the reproduction kernel Q, and analyze its properties. 

The sexual reproduction has two opposite effects on the repartition of the 
population over the phenotypic trait: 

• For each set of two chromosomes, one come from the parent of pheno- 
typic trait t>*, and one from the parent of trait v'^, this process tends 
to give the offspring a trait between v* and v'^. The effect of this phe- 
nomenon is to concentrate the population traits. 

• On the contrary, a variability is maintained in the population by mu- 
tations and recombinations (notice that the effect of recombinations is 
typically much larger than the effect of mutations, see [S]). 

A reasonable assumption is then to assume that in absence of selection, 
a sexual population phenotypic distribution will converge to a given profile 
that only depends on Q. To make this assumption precise, we consider the 
homogeneous sexual reproduction model (J3J) without selection, and with a 
constant population size, that is 

f d t n{t, v) = J J Q(v, v*,v'Jn(t, v*)n(t, <) dv* dv^ - n{t, v) 

\ n(0, v) = n°(v) e L^R), with / n°(v) dv = 1, 1 ' 

then we assume that the long-time dynamics of this model is simple in the 
sense that: 



Assumption 1: 

• for any v*, G R, J Q(v, f*, v'J dv = 1, 
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for any v», G R, J v Q(v,v*,v'^) dv 



2 ' 



• There exists G (from the genetic variance, see [19J) such that for any 
initial population n°, the phenotypic variance of the population con- 
verges to G: 

j ~ J wn °{ w )dw\ n(t,v) dv -» G 2 , as t — > oo, 

• For any initial population n°, the third moment of the phenotypic dis- 
tribution of the population converges to 0: 

J ^v — J wn°(w)dvjj n(t, v) dv — > 0, as t — > oo. 

The reproduction kernel that is typically used in structured population 
models for sexual populations (see [3 [15]) is: 



Q(v,v*,vi) := — ^=e ^ ' • (8) 

We show that Assumption 1 is satisfied for this particular reproduction ker- 
nel: 

Proposition 1 Assumption 1 is satisfied by the reproduction kernel defined 
in ©, with G = V2j 

Remark 1 This assumption is indeed true under the more general assump- 
tion that QiVjVxjV^) = (t * Q(-, v*, v'^jj (v), where f is symetrical with a 

positive variance, Q( v *~^ v * +v, f*, i£) = Q C*^* —v, t>*, v'*) for any v, v#, v'* E 
R, and 

v — Q(v,v*,vjav = C- 



2 I -V \ - 7 - ■» 7 -*/ - j 7 

C < 1. For this and more on this type of problem, we refer to J77| WT^ . 



Proof of Prop. [Q 

The two first conditions of Assumption 1 can be easilly checked. We then 
only prove the two last ones. 
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If we assume w.l.o.g. that J wn°(w) = 0, then 




— fv 2 n(t,v)dv = [ I ( I ' v 2 Q{v,v*,v' if ) dv\ n(t,v'*)n(t,v*) dv*dv[ — f v 2 n(t,v)dv 

>„,.,«-/,.„«. 

= 7 2 — - J v 2 n(t,v)dv. 

We deduce that the variance of n converges to G = for any initial 

condition n°: 

v — w n°(w) dw I n(t, v) dv — >■ G 2 = 2^y 2 , as t — > oo. 



Similarly we have 

fv 3 n(t,v)dv = I I [ I v ? 'Q(v,v*,v'j)dv)n(t,v'j)dVi f dv^— j v 3 n(t,v)dv 





v* + v' # 



n(t, v'„) dv* dvl — / v 3 n(t,v)dv 



-- I v 3 n(t,v)dv, 



4 

so that J v 3 n(t, v) — > as t — > oo. 

□ 

2.3 Rescaling of the structured population model 

To simplify (jSJ), we perform the following rescaling: 



-1/2 

•E @x I T'max ~ I X, 



G 2 


k 


~2V S ~~ 




G 2 


I, 


' 2V S 


2 


G 2 




; 2V S 




V 




= G' 




7 





/ V 1 

max I ' s \ 

~KG \ ' """ ~~ 2V q 1 "■ 



Q(v + * - - ,v*,vi) := GQ(Gv,Gv*,Gvi) , 
where G is given by Assumption 1. Then, fl6]) then becomes: 
d t n(t, x, v) — A x n(t, x, v) 



A A f 

(C — — ) + — (v — Bx) 2 + / n(t, x, w) dw 



n(t, x, v) 




where 



G 2 ( G 2 I 



A . , , I T max 



V, V 2V, 

B : 



G 2 I^ 1/2 



Ga T \" Lax 2K 2 



G 2 J/" 1 



C : 7 [r max 2 ^ 2 j 1. 

Moreover, the rescaled reproduction kernel satisfies the Assumption 1 with 
G = l. 

Remark 2 This is not the only possible rescaling for this equation. As we 
will see in Subsection\2.4\ this particular scaling will allow us to obtain (TT5|) 



as a formal limit of ([9]) when C is large. In fT&J . the authors perform a renor- 
malisation on ([1]) to get only two parameters A and B. Those parameters 
are then defined as: 

A — „, — - , r , B 



2V s {r max - I J 2) ' VW s (r max - I s /2) ' 

The scaling they use wouldn't work here, because they assume that the scaling 
of the variable v that they use does not modify G, the typical phenotypic 
variance of the phenotypic distribution of the population. With (jBJ), G is 
necessarly affected by a scaling in the v variable (see Assumption 1). This is 
also why we couldn't obtain exactly (jTJ), but the slightly different model ( IT3]) 

2.4 Formal limit of the structured population model 

We denote by N, Z, V the following moments of the distribution n(t, x, •): 

N(t,x) :— / n(t,x,v)dv, Z(t,x) :— / v——^ — , -^dv, 



N(t,x) 
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V{t,x) := (v-Z(t,x))' 



,n(t, x, v) 
N(t,x) 



dv. 



We show that iV and Z satisfy the following unclosed equations: 

Proposition 2 Ifn is a solution of fl9]), then the moments of the phenotypic 
distribution of n satisfy: 



d t N(t,x)-A x N{t,x) 

1 + ^(1 - V(t, x)) - ^(Z(t, x) - Bxf - N(t, x) 



N(t,x). (10) 



d t Z(t,x)-A x Z{t,x) 

= 2d x (\og N(t,x))d x Z(t,x) + A{Bx - Z(t,x))V{t,x) 



3 n(t,x,v) 



dv. 



(11) 



Proof of Prop. |2 

To get ffTUj) . we integrate (T9J) along v: 



dtN(t,x)-A x N(t,x) 

d t n(t^v)-AMt^v)d V 

1 + | -^(Z(t,x) - Bxf - N(t,x) 



A 

N(t, x) — — I (v — Bxfn(t, x, v) dv 



1 + |(1 - V(t, x)) - ^(Z(t, x) - Bxf - N{t, x) 



The second equation, (FTP , is obtained as follows: 



N(t,x). 



d t Z{t,x)-A x Z{t,x) 

f n(t,x,v) f n(t,x,v) 

= o t v^rrr. — r- dv — A x I v^rrr. — r- dv 



N(t,x) 



N{t,x) 



N(t,x] 



(d t n(t, x, v) — A x n(t, x, v)) dv 



-(d t N(t,x)-A x N{t,x)) 
+ 2d x {\og N{t,x))d x Z(t,x), 



Z(t,x) 
N(t,x) 
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and then, 



d t Z{t, x) - A x Z{t, x 
v 



N{t,x) 
+ (C + 1) J 



C 



A 



+ —(v — Bx) 2 + / n{t, w) dw ) n(t, x, v) dv 



N(t,x) 




A 



n(t,x,v^)n(t,x,v' ¥ ] 
f n(t, x, w) dw 
A 



1 - -^(Z(t, x) - Bx) 2 -N+—(l- V(t, x))J N(t, x 
+ 2d x {\og N(t,x))d x Z(t,x) 

2d x (log N(t,x))d x Z(t,x) + A(Bx - Z(t,x))V(t,x) 
- ({v-Zf n{t ^ v) 



P *v Q{v, v*,v'*) dv* dv* J dv 



Z(t, x) 
' N(t,x) 



N(t,x) 



dv, 



where we have used the fact that the reproduction kernel does not affect the 
mean phenotypic trait J vQ{v,v', v'*) dv = - ~^ v * . 

□ 

Remark 3 The term 2d x (log N)d x Z is referred to by biologists as the "gene 
flow" term (see f2Mf and [24\ fTPj 26]). This term models the fact that the 
mean phenotype of low density areas are greatly influenced by the phenotypes 
of neighboring high density areas. It is interesting to notice that this term 
does not come from the sexual reproduction term, but from the diffusion term: 



AZ(t,x) = J 



A x n(t,x,v) 



Z{t,x) 



' m+ \ ' dv--^AN(t,x) + 2d x (\ogN(t,x))d x Z(t,x). 
JM{t,x) JM{t,x) 



To close the equations on N and Z obtained in Prop. El notice that 
can be written: 




d t n(t, x, v) — A x n(t, x, v) 

, , n(t,x,v*)n(t,x,v'*) 

+ 
+ 



J n(t, x, w) dw 
— - —(v - Bx) 2 - J n{t, x, w) 
n(t, x, v*)n(t, x, v* 



Q(v, v*,v'*) dv* dv'* — n(t, x, v) 



dw 



n(t, x, v) 




J n(t, x, w) dw 



Q(v,v*,vl) dv* dv* 



so that if C is very large, the first term will dominate the dynamics of the 
population. Since this first term corresponds to the "pure" sexual reproduc- 
tion equation (J7|), and thanks to Assumption 1, it is natural to assume that 
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at all time t > and all locations i£l, 



V(t,x) = J(v-Z{t, x )f^^-dv~G = l, 

Notice that here, thanks to the rescaling performed in the preceding subsec- 
tion, G = 1. 

If we use those properties to close the system of equations on N and Z, 
we get the model: 

f d t N(t, x) - A x N(t, x) = (1 - j(Z(t, x) - Bxf - N(t, x)) N(t, x), 
\ d t Z(t,x) -A x Z(t,x) = 2d x (log N(t,x))d x Z(t,x) + A(Bx - Z(t,x)). 

(13) 

This model is very close to the model ([I]) from [2H US]. The model ([1]) 
was build directly, without the intermediate step of a structured population 
model, and the limits of its applications was unclear (see [22] )■ Our derivation 
shows that the model (fT3"j) is valid (in the sense that it is the formal limit of 
©) if: 

• The reproduction is sexual, 

• The reproduction kernel satisfies Assumption 1, 

• C is large. 

Remark 4 C = 7 (r max — — — 1 is large if the birth rate is large 

compared to the maximal fitness of the population (or many generations are 
necessary to obtain a significant growth of the population, which seems rea- 
sonable in many biological situations) . 

Notice also that for the example of Prop. [21 the convergence of Assump- 
tion 1 are exponentially fast, so that the simplification f |T2l) may be accurate 
even if C is not very large. 

The model ffTBl may hold in other situations, justifying the closure as- 
sumption (fT2|) with other arguments. However, we show in Section H] that 
ffTB"]) cannot hold for asexual populations for the whole range of parameters 
A and B (see Remark [8]). 
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2.5 Derivation of a simplified model 

To simplify (fT3|) . we first apply the following change of variable: 

t := At, x = VAx, 



-Z. 

2 



The rescaled model becomes 



d t N(t, x) - A x N(t, x) = \ (1 - (Z(t, x) - (B/V2)x) 2 - N(t, x)) N(t, x), 
d t Z{t,x) - A x Z(t,x) = 2d x (log N(t,x))d x Z(t,x) + {{B/^2)x - Z{t,x)). 

Now, if we assume that A is very small, N and Z are related by the simple 
relation: 

N(t, x) ~ 1 - (Z(t, x) - (B/V2)x) 2 . (14) 
Therefore we get the simpler model |2|, on Z only: 

O t Z(t,x)-A x Z(t,x) = -4^^) - g/yglgfe^Z 

V ; V ; 1 - (Z{t, x) - [Bj V2)x) 2 v ; 

+ ((fi/v / 2)x-Z(t,x)). 

Remark 5 In fig)/ , t/je range of A that has been considered was A G [0.001, 1]. 
Our approximation assuming that A is small thus seems reasonable. 



Another simplification had been proposed in /TP) /, where the equation on 
N was replaced by 

N := kexp(-f (1- A(Z - Bx) 2 )) . (15) 
dl]) then simplifies considerably: 

d t Z(t, x) - A x Z(t, x) = A{Bx - Z(t, x)) [1 - ^d x Z(t, x)(B - d x Z(t, x))} . 

However, the simplification (1151) seems independent of ([I]). Our simplifica- 
tion has the advantage to rely on a clearer assumption: (|T6|) is the formal 
limit of (|T3l) when A is small. 
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3 Dynamics of sexual populations 



3.1 Well-posedness of the model 

By replacing W = Z - (B/y/2)x in (ED we obtain the following equation 

d T WW r- 
d t W - A X W = -4 * TT/2 (d x W + B/y/2) - W, (16) 
1 — W z 

with —1 < W < 1. This equation has a singularity for W = ±l.The existence 
of singularities is an obstacle to have a well-defined problem. However, as we 
will see in section 13.21 the singularities lead to the existence of propagative 
fronts. In most of the cases in the classical study of propagative fronts, one 
proves the existence of propagative fronts that connect an unstable steady 
state to a stable steady state. Here the situation is different. The propaga- 
tive fronts connect the unstable steady state W = to the singular point 
W — — 1. While the presence of singularities is crucial to observe propagative 
fronts, it is an obstacle to prove uniqueness or comparison results to compare 
the solution with the propagative fronts. Nevertheless we are able to intro- 
duce an approached model where the uniqueness and comparison principles 
hold. 

Since ([TBI) is singular, we approximate it by the following model 
d t W s - A X W S = - \_ w! + s (d*W s + B/V2) - ; _ ^ - - , (17) 
with 

W s (t = 0, •) = <(■)• (18) 

With this choice of approximation we avoid the singularities and transform 
the singularity in —1 into a stable steady state (the stability is for the ODE 
formulation presented in section I3T21 . 
Under the assumption 

- 1 < W° < 1, (19) 

equation (fT7|) has a smooth solution that stays between —1 and 1 by the 
maximum principle. Using the following assumption on the initial data 

\d x W$\ < d, (20) 

with C\ a positive constant, we prove a uniform Lipschitz bound for the Wg's 
and we deduce that the W s, s converge to a viscosity solution of a variant of 
equation f|T6l) (see [TH [2] for general introduction to the theory of viscosity 
solutions): 
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Proposition 3 Under assumptions ( TTTJ|) and f[2"Uj) . we have that the solutions 
of ffTTI) are uniformly bounded and Lipschitz: 

-1 < W s (t,x) < 1, \d x W s (t,x)\ < C 2 , for all (t, x) e R + x R. 

Consequently, after extraction of a subsequence, the Wx 's converge to a con- 
tinuous function W that is a viscosity solution of 

(l-W 2 )d t W-(l-W 2 )A x W = -4d x WW(d x W+B /V2)-(l -W 2 )W. (21) 

We notice that equation f l2Tl) is the original model f|T6|) multiplied by 1 — W 2 . 



Proof of Prop. [3j 

We differentiate equation ffTB]) with respect to x and obtain 

d t d x W s - A x d x W s = -4 T= ^(2d x W s + B/V2)d x (d x W s ) 

7, _l_ R U/o\ 

~a I > 



-id,Wf(%W s + B/V2) ,^tL (22) 



(i - sfflfr) 

where the last term comes from 

f W s (l-W 2 ) \ ( 5W S \ _ ( l-W 2 + 5 + 2W 2 \ 

We multiply fl22|) by d x Ws and devide by |^ x W<j| and obtain 

d t \d x W s \ - A x \d x W 5 \ < -Aj^-^Ws + B/V2)d x (\d x W 6 \) 

-Ad x W 2 (d x W s + B/V2) sgn(d x W s ) 

(23) 

It follows that, for 5 < 1, (c^W^I is a subsolution of the following equation 

d t9 -A x g =a(t,x)d x g+( { ^j^ y2 ) (-4<? 3 + 2^2 Bg 2 + g) , (24) 
with 

Wx r- 

We choose a positive constant C 2 such that — 4Cf + \/2BC 2 + C 2 < and 
Ci < C*2. The constant C 2 is a supersolution of the equation above and thus 

\d x W s \<C 2 . 
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We proved that the Wg's are uniformly Lipschitz continuous in space. More- 
over we know that the Ws's are bounded. It follows that the Ws's are uni- 
formly continuous in time (see [3]). Using Arzela Ascoli Theorem we conclude 
that, after extraction of a subsequence, the Ws's converge to a continuous 
function W. By the stability of viscosity solutions (see [2]), W is a viscosity 
solution of (gID. 

□ 

Remark 6 We can relax assumption (}2"Uj) in Proposition^ This is because 
y = ^= + L is also a supersolution of equation ( 12"%|) . for L = L(B) a large 
constant and 5 small. Therefore we have 

It follows that there is a regularizing effect and the W$ 's become uniformly 
Lipschitz, for allt > t > 0, even if they are not uniformly Lipschitz initially. 

We proved that the equation (j2Tl) has a solution in the viscosity sense. 
Unfortunately the viscosity criteria is not enough to impose uniqueness. We 
give a counter-example below: 

Example. Non-uniqueness for equation f[2Tj) : We have the two following 
solutions to equation (T2TT) : 

Wtit, x) = -1, for all (t, x) eR + x R, 

W 2 {t, x) = -e~\ for all {t, x) G R + x R. 

Here the biological solution is the first one. Because W = — 1 corresponds 
to N = 0. Therefore, if initially W(0,x) = —1 for all x G R, we expect that 
W(t, ■) = —1, for all t G R + . Otherwise some mass is created out of nowhere. 
We can easily verify that if Ws(t = 0, ■) = — 1, we have Ws(t, ■) = —1 for all 
t > 0. Therefore our approximation chooses the biological solution. 

In section 13.21 we study the propagative fronts for this model. To be 
able to compare the solutions with the propagative fronts and to show the 
propagation of the density in space, we need a comparison principle. Unfor- 
tunately as we saw above, the equation (|2"TI) does not have a unique viscosity 
solution and therefore it does not admit a comparison principle. However we 
can prove a comparison principle for the approached model. We first recall 
it's definition: 
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Definition 1 Comparison principle: Equation L(D 2 u, Du,u,x,t) = veri- 
fies a comparison principle, if for any subsolution w 1 and supersolution w 2 
of L such that w 1 (0,x) < w 2 (0,x), we have 

w l (t, x) < w 2 (t, x), for all (t, x) eR + x R. 

We prove that there is a comparison property for (fTTI) . In particular ( TTTj) 
has a unique solution. 

Proposition 4 The problem ( ITTj) admits a comparison principle in the set 
of solutions {-1 < W < 1}. 

Proof of Prop. f4j 

We suppose that W\ and W 2 are respectively subsolution and supersolu- 
tion of ( ITTj) and 

Wi(t = o,o < w 2 (t = o,-). 

We prove that H^i < for all (t, x) G M + x R. Let (t, 2) a maximum point of 
W\ — W 2 . Since it is a maximum point we have d x Wi(t, x) = d x W2(i, x) = p. 
Therefore we have 

d t (W x - W 2 )(i,x) - A(Wi - W 2 )(t,z) 

r- (l + WiWr> + 5) 

* ~ Ap ip + B/ ^ 2) (i-w 2 + m -w 2 + 5) ^ ~ W * & *> 

(1 +5)(1- W 2 - W,W 2 - W 2 ) + W,W 2 + W 2 W 2 _ 



(1 - W 2 + 5){1- W 2 + 5) 

In the previous section we proved that |<?a.VK| is bounded. Thus p(p + B/\r2) 
is bounded. Moreover W% and W 2 are bounded and 

1 - W 2 + 5 > 5, for i = 1, 2. 

Therefore the coefficient of W\ — W 2 is bounded. Following the classical max- 
imum principle we deduce that equation (fTTI) admits a comparison principle. 

□ 



3.2 Existence of propagative fronts and steady popu- 
lations 

We are interested in propagation fronts, that is solutions of (T5]) of the form 
Z(t,x) = (B/y/2)x + U(x — ut). The equation ([2]) becomes: 

-uU' - U" = -i-^L(U' + B/y/2) - U. 
1 — U z 
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Figure 1: On the left, we represent solutions of the ODE defined by the vector 
field (j25|) for B := 0.5 and v = 3.75, 4.75, 5.75 (the vector field represented 
corresponds to v — 4.75). A propagative front exists for those two last values 
of v only, which we represent on the right. 



If we denote by V := U', finding a propagative front is then equivalent to 
find a solution defined on K to the ODE given by the vector field 

Mu,v) = v, 

F V (U, V) = -vV + 4^(V + B/y/2) + U. ^ 
To have a meaning with respect to ffToT) . those solutions must satisfy 

u(t) e[-l,l]. 



Proposition 5 For any B > 0, there exists vb G R such that (jlb|) /ias a 
propagative front of speed v, Z(t, x) = (B/y/2)x + U(x — vt), satisfying 

U(x) as x — >■ — oo, — >• — 1 as x — >■ +oo, 

«/ and on/y if v > vb- 

The propagative front with speed v is unique (up to a translation), and 
vb is a decreasing function of B. 

For each B > 0, there exists a one-parameter family of propagative fronts. 
By analogy to the KPP-Fisher equation, one can guess that there is only one 
stable propagative front, the one with the least speed. Those fronts would 
be invasive fronts if vb > 0, and extinction fronts if vb < 0. The proposition 
[6] shows that steady populations exist in this second case only. 

We notice that, since the model does not admit a comparison princi- 
ple, we cannot use the usual methods used in the study of the KPP-Fisher 



17 



Figure 2: On the left, we represent solutions of the ODE defined by the vector 
field (|25)) for B := 1.5 and v = 0. Two of the three solutions represented 
define steady-states of ffTEj) . which we represent on the right. 



equation, to study the stability of propagative fronts rigorously. Neverthe- 
less, if we choose those solutions of (fl6l) that are obtained as the limits of 
the approached solutions Ws's, and since the appraoched model admits a 
comparison principle, one can expect that the comparison principle be true 
for these limit solutions. The comparison principle would in particular al- 
low us to compare the solutions with the propagative fronts and prove the 
propagation of the population by the minimal spead of propagative fronts. 

Proposition 6 The equation ([2]) has a non-trivial steady-state if and only 
if vb < 0. Steady-states Z(t,x) = (B/\/2)x + U(x) satisfy 

U{x) — )■ 1 as x — oo, U{x) — > — 1 as x — )■ +oo. 

As one can see in the proof, if vb < 0, there exists indeed a whole family of 
steady-states. If we assume that U(0) = (to avoid the translation invariance 
of the problem), then the family of steady-states can be parametrized by 
U'(0) G (-K B ,0), for some K B > 0. 

To show those two propositions, we will use the two following lemma: 

Lemma 1 Let v G 1R. There exist two (up to a shift in the time vari- 
able) solutions (u,v) to the ODE defined by the vector field (1251) such that 
(u, v)(t) — > (0, 0) as t — > — oo. At most one of them is globally defined, which 
satisfies: 



2 l 

For this solution, u is strictly decreasing. 



(u,v)(t) ~ -CL e *^'(l, 
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Lemma 2 Let U\ < v 2 , and (u vi ,v vi ), (u Ul ,v Ul ) be the corresponding solu- 
tions given by LemmaUl If for tx, t 2 £ R, 

/ u vx {t x ) = u U2 (t 2 ) 

\ V Vl (tx)<V V2 (t 2 ), ^ 

then, for any t[ > t x , t' 2 > t 2 such that u Ul (t[) = u U2 (t' 2 ), we have v Ul (t[) < 
v V2 {t' 2 ), and this inequality is strict if Vx < v 2 . 



Proof of Lem. [TJ 

The Differential of the vector field F in (0, 0) is 



DF( 0<0 ) 



1 

1 -v 



and (0,0) is hyperbolic (det DF/q^ < 0). The Hartman-Grobman Theorem 
then applies, and there exists only two (non-trivial) solutions (u, v) satisfy- 
ing (u,v)(t) -^ t -¥-oo (0,0). Since the eigenvector associated to the positive 
eigenvalue of DF( ,o) is (1, v/ ' /2 + 4 ~ iy ) ; those two solutions are equivalent to 



{u,v)(t) ~ t ^-oo ±C±e 2 



for some C_, C< > 0. 



The solution such that (u, v)(t) ~ t _^_ 00 C + e 2 4 ^l, ^ 2 + 4 " \ satisfies 

u(t) > 0, v(t) > for some 1 Since Fu(u(i),V) > for V > v(t) and 
Fy(U, v(t)) > Fy(u(i),v(t)) > for U £ [u{i), 1), the solution cannot escape 
[u(i), 1) x [v(t), 00). In particular, for t > i , u'(t) = v(t) > v{t) > and since 
the vector field is not defined for U = 1, the solution cannot be global. 

The other solution satisfies u{t) < 0, v (t) < for t > i. Since Fu(0, V) < 
for V < and F V (U, 0) < for U < 0, the solution cannot escape M?_, and 
in particular, u'(t) = v(t) < 0, which shows that u is strictly decreasing at 
all times. 

□ 



Proof of Lem. 

We know that u v is strictly decreasing, we can thus define the graph of 
(u Vl ,v Ul ). 
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We assume that t±, t 2 are the smallest points respectively in (t 1 , oo) and 
in (t 2 ,oo) such that (u U2 ,v U2 )(i 2 ) = (u Ul ,v Ul )(ii). We have, 

< <(* 2 ), 

this inequality being strict if z/i < u 2 . It follows that > ^ 2 ^| . The 

graph of (u U2 ,v U2 ) can thus only cross the graph of (m^ , ) from below to 
above, when t increases. This is enough to conclude that v Ul (t[) < v U2 (if 2 ), 
for all t[ > ti and t' 2 > t 2 . Moreover, the latter inequality is strict if v\ < v 2 . 
This completes the proof of Lem. [2j 

□ 



Proof of Prop. [5} 

Step 1: We show that the solution (u,v) given by Lem. [1] satisfies either 
v it) = —B/y/2 for some i, or (u,v)(t) — > (—1,0) as t — > +oo. Moreover, the 
solution is global in this last case only. 

Since u is strictly decreasing and F V (—1,V) = oo for V G (—B/y/2,0), 
only two situations are possible: either v(i) = —B/y/2 for some i < oo, or 
(u, v)(t) (-1, 0) as t t G R U {oo}. 

If t>(t) = -B/y/2, then < -B/y/2 for all i > i. This is because 
F V (U,-B I y/2) = uB/y/2 + U < F v (u(t),-B/V2) < for U G (-l,u(t)) 
and since w is decreasing. Therefore, u'(t) = v(t) < —B/y/2 for all t > t. 
From the latter together with Fy(—l, V) = — oo for V < —B/y/2 we obtain 
that the solution cannot be globally defined. 

Let (u,v)(t) -»• (-1,0) as t ->• £ G R U {oo}. We show that t = oo. For 
([/, V) close to (-1, 0), we have F v (*/, V) ~ "J* ^ - 1. Then, 



d 
dt 



1 + u 



(t) 



F v {u(t), v{t)){l + u(t)) - v{t)Fu(u(t),v(t)) 



1 



1 



(l + «(t)) 2 



> o, 



(H 
-2B 



nit)) 2 
v(t) 



1 (1 + «(*))-„(*)' 



>/2 l + 

ifi v(t) -w(t) 2 - (l + n(t)) 



v/2 
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if v(t) <=^(1 + u(t)) and (u(t),v(t)) is close to (-1, 0). 

Let t be such that (u,v)(t) is close to (—1,0) for t > t. Then, > 

min (ij^j, ^r) for all t > t, and thus, u'(*) = > -C(u(t) + 1), which 
implies the estimate 

u(t)>-l + (u(t) + l)e~ c ^. 

Since (it, — > (—1,0) as t — > t , it follows that t = +oo, and (u,v) is 
indeed globally defined. 

Step 2: We show next that there exists a constant vb such that, there 
exists a propagative front if and only if v > vb- 

If v > we have F V (U, -B/y/2) = - U > for U G (-1, 1). 

Therefore, the solution given by Lem. [1] cannot cross the line V = -B/y/2, 
and thus it defines a propagative front thanks to Step 1. We deduce that, 
there exists a propagative front if v is large enough. 

On [-1/V2,0] x [-B/y/2,% we have, for v < -(1 + 2^2)5, 



V") < -uV - 4V(y + B/y/2) <BV <B F V (U, V). 

It follows that the solution given by Lem. [1] necessarily crosses the line 
V = —B, and thus it does not define a propagative front thanks to Step 1. 
We deduce that the model does not admit a propagative front if — v is large 
enough. 

Consider a solution (u Ul ,v Ul )(t) given by Lem. [T]for some ui, that con- 
verges to (—1, 0) as t — > oo, and v-i > V\. Then, since is a decreasing; 
function of v and 



~ -c_ e ^ — [ 1, y-HL 



the graph of (ti^, v Vl ){t) is below the graph of (u U2 , v U2 )(t) for t << 0. Thanks 
to Lem. 2, this implies that the whole graph of (u Ul ,v Ul ) is below the graph 
of (u V2 ,v V2 ). Using the latter and Step 1 we obtain that (u U2 ,v U2 ) defines a 
propagative front. 

Finally, we show that vb is a decreasing function of B. Firstly we notice 
that, for B 1 < B 2 , we have F// 1 = F//\ and F// 1 < F// 2 on (-1,0] x R_. 
It follows that, thanks to Step 1, if the solution given by Lem. [T] for Bi 
converges to (—1, 0) as t — > oo, so does the one associated to B 2 . This shows 
that vb is a decreasing function of B. 
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□ 



Proof of Prop. 

Assume that vb < 0. Then, for v = ^y, the solution (u UB /2,v„ B /2) given 
by Lem. (TJ is globally defined, and it satisfies {u VB /2,v VB /2){t) — > (—1,0) as 
t — > +00. Moreover we have 



(u UB , 2 ,v VB/2 ){t) ~ t ^_ oc -C_e 2 Ml, v _ 

and > 1 

Consider now the vector field (1251) for v = 0. Since (wo,t>o)(t) 
— C"_e*(l,l), for t small enough, (ito,uo)(i) is strictly above the graph of 
(u UB /2,v UB / 2 ). Let v be such that (u (t),v) is strictly between (u ,v )(t) 
and the graph of (u UB /2,v UB / 2 ). We define (u, u) to be the solution of the 
ODE given by the vector fields (|2"5T) such that (tt, v)(0) = (uo(t),v) and 
^ = 0. Then, thanks to Lem. [2, (u,v) is defined on M + . Moreover, since 
(wo,fo)(^) — >■ (0,0) as t — t- —00 and (0,0) is a hyperbolic point, there exists 
t > such that = 0. 

By symmetry, (u,v)(t + 1) = (—u,v)(t — t), and thus, (u,v) is globally 
defined and satisfies (u,v)(t) —> (±1,0) as t — > ±00. This completes the 
proof of proposition [6j 

□ 



4 The case of asexual populations 
4.1 The model 

We consider here the same fitness ([3]) as we considered for sexual populations. 
If we additionally model mutations through a diffusion of rate (for more 
on the different ways to model mutations, see [ID]), then the evolution of a 
population structured by a phenotypic trait v only can be modeled by the 
classical model (see e.g. [8| fT5]): 



d t n(t,v) - alA v n(t,v) 



1 

r max ~ 



2K 



1 f 

— / n(t, w) dw 
K J 



n(t, v) 



Just as in the sexual case, we add a spatial structure to this model, the 
population is then structured by both a phenotypic trait v as above, but also 
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by a space variable xGK: n(t, x, v). We assume that the selection-mutation 
process described above occurs locally in space, but that individuals move 
randomly in space (which we model by a diffusion of rate a x ), and that the 
trait of optimal adaptation changes linearly in space (see ([5])). We then get 
the following model for asexual populations (-rf represents lethal mutations): 



d t n(t, x, v) — a 2 A x n(t, x, v) — a 2 A v n{t, x, v) 

Is 1 1 f 

= r max - - —{y -bxf - — / n(t,x,w)dw 
We then rescale the problem as follows: 



n(t, x, v). 



—t, n = K(r r 



x 



Trnax 2 



Tmax 2 



and define the two parameters 



.4 



1 



2V s a vX r. 



B := ^b. 



max 2 

Then, we obtain the following rescaled model: 

d t n(t, x, v) — A x n(t, x, v) — A v n(t, x, v) 



1 — A(v — Bx) 2 — I n(t,x,v')dv' 



n(t, x, v) 



(27) 



For an existence theory for this equation, we refer to [H [27]. Notice 
that an integration of (|27|) over the v variable provides the uniform in time 
estimate: 

IMU^^l 1 ^)) — max (II^IU 00 ^ 1 ^)), l) • (28) 

4.2 Qualitative properties of asexual population 

We show that if A(l+B 2 ) > 1, then, for any initial population, the population 
goes extinct when t — > oo: 

Proposition 7 Assume A(l + B 2 ) > 1. For any initial population n° e L°° , 
the population will go extinct exponentially fast when t — > oo: 



iv^-^Vt,,-) 

where c = y/A{\ + B 2 ) - 1. 



0{e 



-ct\ 
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If A(l + B 2 ) < 1, the population survives, and does not remain confined 
in a given part of the space set: 



Proposition 8 Assume A(l + B 2 ) < 1. There exists X, k > depending 
only on A and B such that for any initial population n° ^ 0, any T\ > 
0, Xq G R, there exists T 2 > T\ such that 



Remark 7 We believe that if A(l + B 2 ) < 1, then the population indeed 
invades the whole space in the stronger sense that there exists X, k > 
depending only on A and B such that for any initial population /° ^ 0, and 
any x G R, there exists T > such that for all t >T, 



We were unfortunatly unable to show this stronger result. 

On the contrary, it is not clear that the result of Prop. O would be true 
locally, that is if X could be choosen arbitrarly small. It is indeed known 
that reaction- diffusion equations with integral terms may lead to complicated 
invasion fronts , see f^/. 

Remark 8 Prop. [7| and Prop. show that the dynamics of asexual popula- 
tions is different from the one of sexual populations: For sexual populations, 
numerical simulations from fT9^ show that for some parameters, a popula- 
tion can survive, but remain confined in a restricted area (this observation 
being comforted by Prop. 0, although the model (j2J) is too simple to model 
extinction cases). This shows a property that is well known experimentally: 
asexual populations have a wider geographic distributions than asexual popu- 
lations (see Wh\ Wl). 

Surprisingly, Prop. \^also holds for asexual populations (see also Rem. 
The reason why the model ffT3l) cannot be used to model asexual population is 
that the closure condition fTT2l) is based on the Assumption 1, that is on the 
sexual reproduction kernel, which cannot be satisfied in the case of asexual 
populations. 

Proof of Prop. [TJ 

We consider the following function, for some A > 0: 





<t>(t,x,v) := Ae-tv^^-DV^y^"^ 2 
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that is a strong solution of: 

d t <f>(t, x, v) - A x (f>(t, x, v) - A v 4>(t, x, v)=[l- A(v - Bx) 2 ] <f>(t, x, v). (29) 
Notice that since —n(t,x,v) f n(t,x,v') dv' < 0, n is a subsolution of (T29l) . 

, then n° < 0(0, • 



Moreover if we choose A : - 



to the comparison principle for the equation 



Thanks 



we obtain 



n(t, x, v) < (f)(t, x,v), Vt > 0, x E R, v E 
Since A(l + B 2 ) > 1, <ft vanishes as t tends to infinity: 



< 



< Ce-(V^+ Bi )-Di o, as t -+ oo. 



This concludes the proof of Prop. [3 



□ 



Proof of Prop. [8} 

Step i: We show that \\n(i + 1, -, •)IU'»([xo-x/2,xo+x/2] I iii(*)) can be con- 
trolled by ||n(t, •, OIUmN-x^o+xIxr for some X > 0. 

To show this, we notice that (f(t,x,v) := f-re « is the fundamental 
solution of: 

d t ip{t, x, v) - A x if(t, x, v) - A v (f(t, x, v) = ip(t, x, v)., (30) 

and that n is a subsolution of the same equation. Then, the comparison 
principle for fl30|) shows that 

n(i+l,-,-) < (v(l)* x , v n(t))(-,-), 

since n(i, -, •) < (y?(£) * X}V n°)(-, •). In particular, we have 

\\n(t + 1, -, ^WL^dxo-X^xo+X^lL^iv)) 

< WfO-) *x,v ^(^)||L°°([ :C o-X/2, :E o+X/2],L 1 (f)) 

< 11^(1) *x,v (n(i)t[ X0 -.x,x +X]xM.)\\L°°({x ~X/2,x +X/2],L 1 {v)) 

+ ||y?(l) *x,v (n(t )t[x -X,x +XYxR)\\L^([x -X/2,x +X/2],L^(v)) 

< ||™(* )IU 1 ([ :C o-X, a; o+X]xR) 11^(1, OIU^^.L 1 ^)) 

+ ll"(f)IU°°(a!,L 1 Cy)) f f tp{l,x,v)dvdx 

e 



< 



'\x-x \>X/2 
\n(t)\\Li([x -X,xo+X]xR) 4 



\n(i)\\ 



2e 



L°°(x,L 1 (v)) 



'4:71 '"' " ' ' V47T Jx/2 

<C\\n{t)\\ L 1 ([x -X,x +X]xR) + o x ^oo(l)||n(t )|| 



e 4 dx 
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Step 2: We show that if t > is large enough, then n can be minored by 
a Gaussian function. 
We define: 



ii(t,x,v) : = 



e- ct 



Ant 



-e 44 e 



Then, 

(d t -A x - A v )ijj(t,x,v) 

(a(l + B 2 ) - c) - a 2 (l + B 2 )(v - Bxf 



a- 



(v - Bxf 



It 



ip(t,x,v) 



< [(«(! + B 2 )-c) - a 2 (l + B 2 )(v - Bxf]ij(t,x,v). 



Let xq G R, and R > large enough for J n° R (x,v)dxdv > to hold, 
where 



n 



R - ^{(z^lz-rroP + M 2 ^.?? 2 }- 



We define n(t, x, v) := (ip(t, •, •) * Xj „ n^) (x, t>). Then, 

(d t - A x - A v )n(t,x,v) 

= ((d t ^(t, •) - A x ^(t, •, •) - A v 7j;(t, •)) **,„ <)(t, x, u) 

< ([((a(l + 5 2 ) - c) - a 2 (l + 5 2 )(t; - Bxf) ^(t, •, •)] *«,„ n° ) (t,x, v ) 

< (a(l + £? 2 ) - c) ra(*, u) - a 2 (l + 5 2 ) (((u - Bxf^j(t, •, •)) <) (t, x, u). 

we can estimate the last term using the fact that supp n° R is bounded, and a 
Young inequality: 

- (((v - Bxfip(t, •)) * x ,v n° R ) (t, x, v) 

= - J J ((v - Bx) - (V - Bx')f ip(t, x-x',v- v')n° R (x', v') dv'dx' 

= - J J ((v- Bxf + (v' - Bx'f - 2(v - Bx){v' - Bx')) 

ip(t, x — x',v — v')n R (x', v') dv'dx' 

— [~( v ~~ Bx f + \v — Bx\] n(t, x, v) 

I 



< 



-(1 - 5)(v - Bxf + - 



n(t, x, v) 



Finally, we get that: 

(dt - A x - A v )n(t,x,v 
< 



a(l + B 2 )(l + -) - c - (1 - 5)a 2 (l + B 2 )(v - Bxf n(t, x, v), 
o J 
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We now choose a and c as follows: 



a 



A 

(1 + B 2 )(1-S) 



c = a(l + B 2 )(l + j) +max(||n°|| L <»( a , > jr,i( t) )), l) 



Ail + B 2 ) A 
\_ 6 (1 + B 2 ) + + max (IKIUo^,^)), 1) , 

then n is a subsolution of: 



<9t^>(£, x, u) - A x if)(t, x, v) - A v 4>(t, x, v) 

= [0 - A(v - Bxf - max (lln ^^^)), l)] ip(t,x,v). (31) 

Since IMIl^^^l 1 ^)) < m a x (H^IU^^i 1 ^))) l), n is a supersolution of (jHTl) . 
and since n° > the comparison principle for (13"T1) shows that n > n, for 
all £ > 0, and in particular, 

„(*,„,)> -»( I - I ,) tI , t .)/ „&(*,, +W ,«0*«to 

ll(l/>*o)ll<-R J\\(y,w)\\<R 

e C * / (a:-3: + y) 2 + ( t »+ t u) 2 (u + m-BQ -xg+y) 

> C mm e 4 * e 2 



4vri ||(y,w)||<Ji 

> exp | - ^-((x - x ) 2 + R 2 + 2R\x - x \ + (v - Bx) 2 

+2(5|x - x | + £|x | + R)\v - Bx\ + B 2 (x - x ) 2 
+2B(B\x \ + R)\x - x \ + (B\x \ + R) 2 )} 
exp {-a ((u - 5x) 2 + 2(R + BR + B\x \) \v - Bx\ + (R + BR + B\x \) 2 ) } 

and then, thanks to Young inequalities, we get, for any 8' > 0: 

e -ct r Q 
n(t,x,v) > C s , — exp i - — (x - x ) 2 

exp |- + 8' + (v - 5x)' 

where depends on 5', Xo, B, f^ y w ^\ <R n °( x o + V, W ) > 0. Finally, for any 
xq G R, /ii, //2 > and Ti > 0, there exists r > Ti, A = A(r) > such that 



n(r, x, v) > Xe nVi+B* w 7 l ' e ^ 2 



Vx, d G E x 
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Step 3: We show the result. 

Since A(l + B 2 ) < 1, there exists \i > such that (this defines the \i and 
then the k): 

1 - yjA{\ + B 2 ) - (i(3 + B 2 ) > 0. (32) 

We define: 

X (t,x,v) := e (i-^^-M3 +B 2 ))t e -K^+^)(-^) 2 e -^ 2 ) 

which satisfies: 

<9 t x(£, u ) - & x x(t, x, v) - A,,x(*, x, v) 

1- fi - (VJ + fiy/l + B 2 ) 2 {v - Bx) 2 - /iV x(*,a;,u) 

< [l - fi - A(v - Bx) 2 - max (H^Hl^l^))' X ) l \x-*o\>x] x(t,x,v), 
if X is large enough. Then, x is a subsolution of (13"3"|) : 

<9 t x(t, £, v) - Aa.x(*, a?, v) - A„x(t, a:, «) 

= [l- ii- A(x - Bv) 2 - max (||n || it x> (a . )L i (v)) , l) l\ x - Xo \ >x \ x(t,x,v). 

(33) 

Thanks to Step 1, provided that X is large, there exists k > such that if 
1 ([x -x,x +x]xR) — K for all times t > T 1; then, for any i > Ti + 1, 
n is bounded as follows: 

n(t, -,v) dv < ji. 

L°°([x -X/2,x +X/2]) 

Then, 1 1-> n(7i + 1 + t, •, •) is a super-solution of fl3"3"|) for £ > 0. 
Thanks to step 2, there exists r > 7\ and A > such that 

n(r, x, v) > A x(0, x, t>), Vx, v e K x E. 

Then, thanks to the comparison principle, n(r + i, •, •) > \x(t, •, •) for t > 0, 
and in particular, 

||n(T + t, v )IU> > A||x(t, v )IUi 



thanks to ( 13"2"j) . This shows that the assertion that \\n(t, •, OIUHN-x^o+xixi 
k for all times t>T\ cannot be valid, and shows Prop. [SJ 



< 



□ 
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